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ABSTRACT 


The forecast method introduced by 


Korsos et al. 


fj2014j) is generalised from 


the horizontal magnetic gradient {Gm), dehned between two opposite polarity 
spots, to all spots within an appropriately dehned region close to the magnetic 
neutral line of an active region. This novel approach is not limited to searching 
for the largest Gm of two single spots as in previous methods. Instead, the pre¬ 
hare conditions of the evolution of spot groups is captured by the introduction of 
the weighted horizontal magnetic gradient, or WGm- This new proxy enables the 
potential of forecasting hares stronger than M5. The improved capability includes 
(i) the prediction of hare onset time and (ii) an assessment whether a hare is 
followed by another event within about 18 hours. The prediction of onset time is 
found to be more accurate here. A linear relationship is established between the 
duration of converging motion and the time elapsed from the moment of closest 
position to that of the hare onset of opposite polarity spot groups. The other 
promising relationship is between the maximum of the WGm prior to haring and 
the value of WGm at the moment of the initial hare onset in the case of multiple 
haring. We found that when the WGm decreases by about 54%, then there is 
no second hare. If, however, when the WGm decreases less than 42%, then there 
will be likely a follow-up hare stronger than M5. This new capability may be 
useful for an automated hare prediction tool. 


Subject headings: Sun: flares — sunspots 
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Introduction 


The endeavour to establish reliable flare forecast methods has resulted in numerous 


attempts to identify promising physical quantit ies, proxies, 


patterns to predict an imminent flare (see e.g. 


Ben 2 


Sawyer 


1986 


eatures and behavioural 


Hochedez et al. 


2005 


20081 ). The ultimate task is to construct a diagnostic tool that determines the 


unique conditions le ading to flaring events, triggered by magnetic reconnection (see e.g. 


Yamada et al. 


20101) . The empirical study presented here focuses on features at the solar 


surface, namely investigating the pre-flare dynamics of sunspot groups. 


Most of the attempts that developed flare forecast tools employ suitably defined 
(and derived) quantities from magnetograms. The majority of previous efforts studied 
the behaviour of the horizontal gradient of the line-of-sight component of the magnetic 
held, usually a well-observed property of an active region (AR). By investigating the Solar 
and Heliospheri c Observatory’s Michelson Doppler Imager (SOHO/MDI) magnetograms, 
Schriiverl (120071) found that energetic flares are connected to the separation line of opposite 
pola rity regions with high mag netic gradient. Further works followed suit, including 


e.g. 


Mason &: Hoeksemal (120101) who studied the gradient-weighted inversion line length 


which exhibits a significant increase prior to flares. Th e maximum hori z ontal gradie nt 


and t h e length of t 


(l2006h. 


r e neu tral l ine were c onside red by 


Huang et al 


Criscuoli et al. 


(l2010[ ) and 


Yu et al. 


( 20091) . however. 


Cui et al. 

(2006 

), 

Jing et al. 


(I 2 OIOI). The fractal structure was addressed by 


GeorgoulisI (120121) was sceptical about this approach. In 


the literature, there m e a number of other indicators proposed to measure non-potentiality. 


Leka fc BarnesI (120071) suggested 8 categories of different distributions with an impressive 


list of 29 types of variables defined on them. All these methods above (and others) 
have advantages and caveats, however, none are yet accepted as a universal and reliable 
prediction tool. 
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Surprisingly, sunspots themselves are rarely considered as possible holders of 
information for flare forecast. Arguably, the simplest way to ap ply sunspot data is to 


approach non-potentiality by using their McIntosh classihcation flBloomheld et al. 


2012 ). 


However, the McIntosh classes are only based on morphology, and they do not contain 
perhaps the most important information relevant to th e present context, i. e. the distribution 


of opposite magnetic polarities. Our preceding study (IKorsds et al. 


2014 hereafter Paper 


I) is the hrst attempt to track the details of the evolution of sunspots (umbrae) in order to 
identify signatures of flare imminence. Paper I implemented a proxy horizontal magnetic 
gradient, Gm, dehned between two single spots of opposite polarities close to the magnetic 
neutral line. The pre-flare behaviour of this proxy quantity exhibited characteristic and 
unique patterns: steep rise, high maximum and a gradual decrease prior to flaring. These 
properties may yield a tool for the assessment of flare probability and intensity within a 
2-10 hour window. 


The method developed here allows us to elaborate on some of the most important 
properties of an imminent flare; its intensity and the onset time. The prediction of 
intensity is found to be more reliable, as a linear relationship has been found between the 
pre-flare maximum of Gm and the peak intensity emitted in the 1-8 A range, according to 
Geostationary Operational Environmental Satellite (GOES) x-ray measurement j^. This 
result may be considered to be an indicator of the existing relationship between the proxies 
of free and released energies. Next, the onset time prediction is found to be somewhat less 
precise, as the most probable time of flare onset is between 2 and 10 hours after the Gm 
maximum. The aim of the present work is to hnd more reliable forecasting methods for the 
onset time as well as the likelihood of consecutive flaring. 


^http://www.ngdc.noaa.gov/stp/satellite/goes/dataaccess.html 
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Method of Examination 


)lf|i 


The empirical basis of this study, similar to Paper I, is the SDEo (SOHO/MDI-Debrecen 
Data) sunspot catalogue, the most detailed of its kind covering the years of MDI operations 
(1996-2010). In addition, the intensity peaks of the examined flares the GOES solar flare 
databasqfl is employed. The SDD is efficient for tracking the internal dynamics of spot 


groups. It contains c 


1.5 hrs (iGvdri et al. 


at a o n position, area and magnetic field for all spots with a cadence of 


mil). 


Let us now consider the involvement of all magnetic spots in an appropriately selected 
area at the region of a Polarity Inversion Line (PIL). When a new spot (min. 3 MSH) 
emerges close (within 40 ±5 Mm, which is always enough to capture the emergence) to 
existing spots of opposite polarity, we determine the maximum Gm between the emerging 
and existing spots. Once a pair of spots of opposite polarities with a maximum Gm is 
found, we compute the area-weighted centre between them. Next, there is a defined circular 
area, around this weighted location where Gm is highest, whose diameter is 3°±0.5° in 
Carrington heliographic coordinates. The center of the circle is fixed and spot groups within 
this area are now monitored. 


We assume that the underlying process driving a flare is a collective one between 
nearby spots. Therefore, a new proxy parameter, a generalisation of the Gm, called the 
weighted horizontal magnetic gradient (WGm) is defined to account for this collective 
behaviour: 


WGm = 


E 


R ■■ A ■ 


- E, Bnj ■ A„j 


^pn 


^http://fenyi. solarobs.unideb.hu/SDD/SDD.html 
^ http://www.ngdc.noaa.gov/stp/satellite/goes/dataaccess.htm 


( 1 ) 
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Here, B (determined by f{A) in Paper I) and A denote the mean magnetic field and 
area of umbra. The indices p and n denote positive and negative polarities, i and j are 
their running indices in the selected spot cluster and dpn is the distance between the 
area-weighted centers of two subgroups of opposite polarities in this cluster. 

We have analysed 45 single and 16 multiple flare cases which are stronger than M5 
between 1996 and 2010. This limitation is based on the findings of Paper I, i.e. the present 
method seems to be suitable for energetic flares above M5. Besides some similarities 
between the behaviour of Gm and VPG'm, the current approach results in significant, 
previously unseen pre-flare behavioural patterns. Let us demonstrate the key features with 
two randomly selected, but typical, examples. 

Figure [T] shows a typical active region, AR 8771, with a single flare. The right-hand 
panels of Fig. [T] are: the white light image (top), magnetogram (bottom) and a cartoon 
reconstructing the AR from the SDD catalogue (middle). The WGm (top left) shows a 
steep rise and high maximum (called WG^°'^) followed by its decrease until the flare. 
However, and most importantly, we found characteristic and appealing differences from 
the result of the single spot-pair method, as demonstrated by e.g. the distance diagnostics 
panel (middle left) of the spot groups of opposite fluxes. This plot contains a conspicuous 
dip, indicating a duration of converging-diverging motion of the area-weighted centres that 
seem to be indicative of the next flare for all cases we investigated. The total unsigned 
magnetic flux (bottom left) in this part of the AR shows some increase before the flare but 
it has no special, identifiable characteristic feature. Additionally, we tested a large sample 
of non-flaring spot groups with opposite polarities to determine whether this behaviour 
is found in these cases as well. We found no such behaviour. Therefore, the evolution of 
the distance between the area-weighted centers of spot groups of opposite polarities shows 
potentials for flare forecast. 


Figure [2] is another typical example, i.e. AR 9393 also examined in Paper I, but now 
with multiple flare activities linked only to Area 2 because the highest variability of WGm 
is in Area 2, while Area 1 and 3 do not show such a property (see in Paper I). A series of 
flares is considered multiple, if after the occurring flare there is another flare, within an 
18-hr window, belonging to the same preceding rising phase of WGm and at the same time 
belonging to the same decreasing phase after WG^°‘^ (see e.g. the XI.4 flare followed by 
the X20 flare in the upper left panel of Fig. [2]). The behaviour of WGm is analogous to 
that of the single spot-pair (compare to Gm of Fig. 2 in Paper I), showing a steep rise, high 
maximum and a slower decrease prior to the XI.7 single flare (29 March), and, in spite of 
the limited temporal resolution, arguably a similar pattern is found leading to the follow-up 
ndependent multiple flares (i.e., XI.4 and X20, 2 April). Next, the distance diagnostics 
(middle left) panel of the spot groups of opposite fluxes do show similar conspicuous 
dips as in Fig. [1] Note, these multiple flares are likely unrelated to the single flare as 
they do not fall within the 18-hour window we propose as a requirement for flares to be 
connected. Again, this pre-flare behaviour pattern was not identifiable in the distance 
diagram of the single spot-pair method in Paper I (see its Fig. 2, left column, middle 
panel). Finally, we could not conclude any special and immediate pre-flare behaviour on a 
daily time-scale by investigating the variation of unsigned flux within Area 2 (bottom left) 
because it was decreasing before the XI.7 flare while it was increasing before the XI.4 flare. 
A statistical study also shows that the temporal variation of the unsigned flux does not 


ex. 


of 


libit characteris tic pre-flare signatures. This, however, does not contradict the findings 


Schrijverl (120071 ) about a statistical relationship between the likelihood of X- or M-flares 


and the unsigned flux at SPIL (Strong-gradient Polarity Inversion Line) because selected 
actual states of ARs were investigated and not the time profile of the examined quantities. 


By the above examples of case studies, we are now encouraged to conduct a statistical 
study to confirm (or refute) the above findings for all cases of flaring spots with a strength 










of over M5 available in the SDD catalogue. 


3. Diagnostic potentials with spot dynamics 

We test the proposed diagnostics on a statistical sample by applying the following 
requirements: Firstly, the examined pre-flare variation is within ±70° from the central 
meridian to avoid geometrical foreshortening close to the limb. Next, the flare onset is no 
further eastward from the central meridian than —40° to have sufficient time to follow the 
development of WGm- 

In Paper I, a linear relationship was found between the maximum value of Gm 
preceding a flare and the peak intensity of flares. This behaviour is conhrmed for WGm 
by Fig. [3l 45 single flares (crosses, left panel) and 45 single with additional 16 largest of 
multiple flares (circles, right panel) show a linear relationship between WG'^°‘^ and the 
corresponding GOES flare intensity. Here, we restrict the empirical analysis for flares 
between M5 and X4 classes only. 

Next, the new method revealed further important connections, which is conspicuous, 
in Figs. [Hand[21 This connection is between the durations of converging-diverging motion 
of the centers of opposite polarities. This intriguing pattern was found in all 61 cases 
investigated here. The question rises, whether there is a relationship between the duration 
of the converging motion (the duration from the moment of the hrst point when the distance 
began decreasing to the moment of the minimum point of the parabolic curve) and the 
time elapsed from the moment of minimum distance until the flare onset (duration of the 
diverging motion and the follow-up time until the flare onset). To determine these two time 
intervals for each flare, parabolic curves were htted to their distance data. For a sample see 
the top left panel of Fig. [Tj which is a parabolic £t to the distance data from the left middle 
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panel of Fig. [H showing the converging-diverging behaviour of this relative motion. 


Figure m gives further insight into the relation between these intervals by plotting the 
time from the moment of minimum distance to the flare onset as a function of the duration 
of converging motion. First, the upper right diagram depicts the duration of diverging 
motion as function of the duration of converging motion for the 45 single (crosses) and 16 
multiple (circles) flare cases. Note that the duration of diverging motion is shorter than 
the time period from the moment of minimum distance to that of the flare onset (see 
bottom panels of Fig. |4]). How ever, the converging-m otion phase and the diverging-mot ion 


phase have the same duration. 


Yamada et al. 


(120lOl) found similar properties in laboratory 


reconnection experiments, and called it the ’’push and pull-mode”. The present observations 
are a conhrmation of the laboratory experiment. 


The lower diagrams plot the time from the moment of closest position to the flare 
onset as function of the duration of the converging motion phase. In the cases of multiple 
flares, we investigated the time from the moment of closest position to the first flare onset 
as function of the duration of the converging motion phase. The left/right panel contains 
those cases when the spot groups are younger/older than three days at the time of flare 
onset. The regression lines of the two cases are, surprisingly, different. By estimating the 
time the magnetic helds younger than three days should be distinguished from the older 
ones, the relevant formulae are given in the lower panels of Figure 01 One may be able to 
estimate a rough onset time of the flare. Note the considerable dispersion. If the study area 
is younger than three days then about a mere hour is needed to be added to the duration 
of the corresponding scaled duration of converging motion, where the scale-factor is 1.3 for 
younger ones, to obtain the flare onset time. However, if the area is older than three days 
then the scale-factor is 0.85, and one needs to add 12 ±3 hours to the scaled duration of 
converging motion. 
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Next, a relationship is found between the values of WGm at its maximum prior to 
flaring and at the time of flare onset as visualised in Figured We 

investigate separately the 45 cases when only a single flare took place after the maximum of 
WGm, and the 16 when more flares erupted after the maximum within an 18-hour window 
after the occurring flare on the decreasing phase of the WGm- The left panel depicts the 
cases of single energetic flares (crosses). The right panel depicts the first flares (circles) of 
those ARs where multiple flares are produced. The plots are interpreted as follows: a single 
flare erupts when the WG^^"^ was less than 5 • 10® Wb/m and the WGm decreases by more 
than half of the WG"^^, in the pre-flare phase (for example the XI.4 flare in the NOAA 
8771 on Fig. [T]and the XI.7 flare in the AR NOAA 9393 on Fig. [2]). This is likely due to 
the fact that the magnetic energy in the region decreases so signihcantly during this hrst 
flare that there is simply not enough energy left to release another flare. In that case, if the 
decrease is smaller than half (about 42%) by the onset of the hrst hare, some further haring 
may be expected, meaning that the disturbance of this hrst hare could be forcing opposite 
polarity helds together in the solar atmosphere leading, for example to the ‘homologous’ 
haring that is often observed. On the other hand, the WG"^^ larger than 5 • 10® Wb/m 
seems to be enough in itself to predict a multiple haring event. We cannot comment yet on 
further hares (i.e. third, etc.) in the case of multiple hares as the temporal resolution of the 
SDD catalogue is too coarse. For an example, see the case of AR NOAA 9393 where an X 
1.4 hare was followed by an X 20 one within 12 hours (Fig. [2]). 


4. Discussion 

In this paper, we present advancements in the classihcation of pre-hare conditions with 
an application to hare prediction. 61 cases were investigated in the vicinity of PILs of ARs. 
We assumed the hare onset to be a collective response to their dynamics. This assumption 
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needs further investigations both observationally (with higher resolution) and theoretically 
(e.g. numerical simulations). 

First, we found that the pre-flare behaviour of the weighted horizontal magnetic 
gradient {WGm) exhibits similar patterns to those found with the single spot-pair method: 
steep rise, high maximum and gradual decrease until the flare onset (Figs. [I]12]). Next, 
Figure [3] corroborates the relationship, reported in Paper I, between the maximum of 
Gm and the intensity of single flares. Here, this relationship is modelled as linear one, 
however, the dispersion is considerable and theoretical (e.g. numerical) modelling may 
be necessary to conhrm or refute this relation. There may be a yet unknown physical 
parameter, therefore not accounted for, that would reduce the dispersion. This is the reason 
of restriction on the currently considered GOES classes. However, the relationship found 
can be still regarded to be a link between the proxy measures of the free energy and the 
released energy. A shortfall of the single spot-pair method was that one could not deduce 
the flare intensity in the case of multiple flares. This is now rectihed by the introduction 
of the weighted horizontal magnetic gradient. The spot-group method is now capable of 
providing a rough estimate of the expected largest flare intensity from WG'^^ (Fig. [3]). 

In addition, this method also gives a better estimate of the expected time of flare onset 
(Fig. 0]), and, may be able to predict whether an energetic flare after WG'^^ is the only 
flare or further flare events can be expected (Fig. [5]). 

Let us assume that WGm is a proxy of the available non-potential (i.e. free) energy to 
be released in a spot group. In this case, we may conclude from Fig. O (i) if the maximum 
of the released energy may be over half of the maximum of the accumulated (free) energy, 
no further energetic flare(s) can be expected; (ii) If the maximum of the released flare 
energy is less than about ~42%, further flares are more probable. In short. Fig. fallows 
us to track the variation of the energy balance of ARs and to assess the probabilities of 


12 


consecutive flares and their intensities. 


Last but not least we provide some notes on the estimate of the onset time of an 
imminent flare. Here, its determination is rehned. Paper I only presented a statistics that 
60% of observed energetic flares are between 2-10 hrs after the maximum of Gm- Figure HI 
however, allows a much stronger statement on the expected time of onset due to WGm- 
The hgure uncovers the relationship between the duration of the converging motion of 
opposite polarities (their compression) and the time elapsed betw een the closest posit ion 


and flare onset following the diverging motion (see earlier motions 


Yamada et al. 


2010|). By 


determining the duration of the converging motion the flare onset can now be assessed for 
all cases. We also found that the data points of the motions of younger spot groups have 
smaller dispersion (left of Fig. |4]). 
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Fig. 1.— NOAA AR 8771, for Nov. 23-26, 1999. Right column: continuum white-light image 
(top), reconstruction from SDD (middle), magnetogram (bottom). Left column: variation of 
WGm (top), distance between the area-weighted centers of the spots of opposite polarities 
(middle) and unsigned flux of all spots in the encircled area (bottom). 
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Fig. 2.— Same as Fig. [T]but of NOAA AR 9393 with a single (i.e., XI.7) and multiple (i.e., 
XI.4 and X20) flares, for Mar 26 - 03 Apr 2001. 
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Fig. 3.— GOES intensity of flares as function of the maximum WGm- The left diagram 
shows the intensity of single solar flares (crosses) that occurred after the maximum WGm- 
The right diagram depicts, in addition, the flares (circles) with largest intensity within an 
18-hr interval after reaching WG^^. 
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Fig. 4.— Upper left: The converging-diverging pattern with a parabolic ht to the data for 
AR 8771. Upper right: Relationship between the durations of converging and diverging 
motion for all 61 flares. Crosses (45)/circles (16) indicate single/hrst of multiple flares. 
Please note that apparent number of point may not be 45/16 in the plots because there are 
overlapping data points. Lower panels: Duration of diverging motion until flare onset as 
function of duration of the compressing phase of motion of opposite polarities. The study 
area is younger (left)/older (right) than three days. 
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Fig. 5.— The weighted horizontal magnetic gradient at flare onset as fnnction 

of the maximnm of the weighted horizontal magnetic gradient prior to flare Left 

panel; cases of a single flare. Right panel: first events of mnltiple flares after WG'^^. 




